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NORMAL MODE ANALYSIS OF A ROTATING GROUP OF 
LASHED TURBINE BLADES Y SUBSTRUCTURES 
By A. W. Filstrup 

Westinghouse Research and Development Center 


SUMMARY 


A group of 5 lashed Identical steam turbine blades is studied through the 
use of single level substructuring itsing NASTRAN Level 15.1. An altered 
version, similar to DMA? Program Number 3 of the NASTRAN Newsletter, of Rigid 
Format 13.0 was used. Steady-state displacements and stresses due to 
centrifugal loads are obtained both without and with consideration of 
differential stiffness. The nonrial mode calculations were performed for 
blades at rest and at operating speed. Substructuring lowered the computation 
costs of the analysis by a factor of four. 


INTRODUCTION 


Triangular plate elements have been used by Westinghouse and others (see 
Ref. 1) in NASTRAN to analyze rotating turbine blades. 

There was a need to analyze a group of five lashed 0.79-m (31-lnch) steam 
turbine blades for operation at 60 revolutions per second. Steady-state 
displacements and stresses were needed as well as the natural frequencies, mode 
shapes, and stress patterns. 

Based on NASTRAN calculations on a single 0.79-m blade with associated 
lashing wires, it was decided that a finite element mesh of 700 CTRIA2 elements 
and 407 grid points would be usud to represent each turbine blade. The root 
flexibility was approximated by 11 CELAS2 elements. 

It was discovered that /ipproximately two hundred degrees of freedom would 
be required in the a-set for each blade using Guyan reduction, if accurate 
stress results were to be found for t'ne modes to be evaluated. Whether or not 
Guyan reduction was to be used and whether the inverse power or Givens method 
were used for the eigenvalue extraction, it was apparent that calculation costs 
would have been prohibitive if substructuring were not used. 

This pap It: describes the successful substructuring analysis of the group 
of blades. Th» steady-state stresses were obtained for operation at 60 
revolutions per second and the natural frequencies were obtained for the first 
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nine inodes at both 0 revolutions per second and at 60 revolutions per second. 


METHOD OF ANALYSIS 


The finite element mesh used to represent a turbine blade or substructure 
Is shown In Fi^. 1. The middle sections of the lashing wires and airfoil are si 
Each lashing wire actually resembles u variable thickness plate more than a win 
The auxiliary program which produces these plots views normal to the middle 
surface of the lashing wire. The airfoil Is highly twisted, and near the base 
it Is highly curved. No one viewing angle could provide a clear representation 
of element layout. Thus the auxiliary program which produces the element 
geon^try and isostress plots opens up each cross section. Different scales are 
used for the lashing wires than for the airfoil in figure 1. 

Each blade or substructure has 407 grid points. The 2442 degrees of 
freedom associated with these points are reduced through single point constraint' 
and omits to an a-set of 301 degrees of freedom. One hundred twenty of rhese 
degrees of freedom are at the tips of the lashing wires and are required for 
connecting adjacent blades or substructures. Sixty degrees of freedom are 
common between adjacent blades. 

The combined matrices for the group of five blades then has 1505 less four 
times sixty or 1265 degrees of freedom. Single point constraints to remove 
rotations about the normals to the surface of the exterior lashing wire's reduce 
the system of equations to be solved to 1245. The half~bandwldth Is 301 with 
no active columns. No secondary Guy an reduction was performed to reduce the 
number of degrees of freedom as the resulting bandwidth would have to be 
significantly larger than 301 for accurate results. The inverse power method 
with shifts was used to solve the eigenvalue problems. 

The identical substructure concept as described in Sec. 1.10.5 of refer- 
ence 2 was used. Five phases were required as shown In figure 2. In some Cases 
It was deemed advisable to use more than the one user tape shown between phases. - 
Even though the differential stiffness would be somewhat different for each of 
the five substructures, only one (the center blade) was generated In Phase III 
and used In Phase IV. This approach reduced the total calculation costs by 
about 20%. The i.»mccuracles of this approximation were fslt to be about the 
same as those due to scmie of the other approximations made. Runs IV and V were 
split Into several parts to enable shorter individual runs. 

The mesh for the airfoil was generated by a preprocessor computer program. 
The meshes for the platform and the two lashing wires were generated by hand. 

The Isostreas llnea for the centrifugal loading and for the scaled eigenvectors 
for the airfoil and lashing wires were plotted with a postprocessor computer 
program which reads Images of punched element stress cards. A 

STRESS (PRINT .PUNCH) » ALL 

card was placed In the Case Control Deck. However, job control cards were used 
to store the card Images on t%ro disks and to prevent the punching of cards. 

Over two hundred thousand card Images were produced In the Phase V runs. An 
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Intermediate program was written to enable the isostress plotting program to 
handle the stress information on the disk, more efficiently. 

Table 1 shows the calculation times for the svibstructuring analysis for 
each phase. The mesh was generated on a CDC 6600 computer and the other runs 
were made on an IBM 370-165 computar. 

Table 2 shows the projected calculation times for the analysis of five 
blades without substructuring provided enough disk space were available which 
is extre.:iely doubtful. In addition, checkpointing and restarting would be 
essential due to the extremely long total running times. However, Level 15.1 
NASTRAN requires that this be done on a single physical tape which obviously 
would not hold enough information. Tlie user would be required to use DMAP 
statements to transfer data from one run to ;inother on user tapes rather than 
checkpoint tapes. Even then, some matrices might be too large to fit on a 
single tape. 

When costs of the CALCOMP plotter are added to the computer costs shown, 
the total cost for a nonsubstructuring analysis, if possible, would have been 
about four times the total cost of the substructuring analysis performed in 
this study. 

The arrangement of the NA!>TRAM decks including the Executive Control Decks 
are shown in the appendix. 


RESULTS AND DISCUSSION 


The natural irequencles, mode shapes, and stresses for the first nine modes 
of a group of five lashed rotating steam turbine blades were f *md. The 
natural frequencies, in general, agreed well with experimental values. 

A Campbell Diagram was p^^epared to determine possible resonances during 
various operating conditions. 

The pseudo steadystate deformations and stresses due to the centrifugal 
forces at operating speed were found. This enables the calculations <•! the 
fluid flow through the ro»> of blades through the passages that actually occur 
in operation and not through the passages in the undeformed condition. Thus, 
NASTRAN provides the designer of flexible turbine blades with a tool to help 
obtain near optimal fluid flow characteristics between the airfoils. 

A sample Isostress plot for one of Che surfaces for one of the blades for 
one of the modes is shown in figure 
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RECOMMENDATIONS 
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1. NASTRAN Level 15 with its substructuring capability can and should be 
used for many structural problems. 

2. When preparing data for large problems, a mesh generator computer 
program should be used as much as possible. 

3. For very rigid rotating turbine blades or blade groups. Rigid Formats 1 
and 3 will give accurate results and should be used. For more flexible blades. 
Rigid Formats 4 and 13, which include the differential stiffness matrix should 
be used. For even more flexible blades, it may be necessary to ALTER the 
centripetal acceleration matrix (see Ref. 3) into Rigid Formats 4 and 13. 

4. In order to encourage more users to use the substructure capability of 
NASTRAN and in order to reduce the effort of the user in creating and checking 
DMAP packages and substructuring data, it is urged that substructuring be made 
more automatic (see Ref. 4) . 

5. Rigid Format 13 should be documented in the NASTRAN documentation. 
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CONCLUSIONS 


1. The determination of the natural frequencies, mode shapes and states of 
stress for lashed rotating and non^rotating steam turbine blades is feasible 
using the general purpose computer program NASTRAN. 

2. Substructuring can greatly reduce the computer costs of large problems. 
For the analysis performed here, the total computer expenses including mesh 
generation and stress plotting were one-fourth what they would have been 
without substructuring. The NASTRAN runs cost one-sixth as much using 
substructuring than they would have cost without substructuring. 

3. Choice of the proper root flexibility is important to produce accurate 
frequencies and stresses for all modes. 

4. Mode shapes and isostress lines for the fifth through ninth modes 
varied significantly between those found at 0 revolutions per second and those 
found at 60 revolutions per second. This variation is due both to the flexi- 
bility of the blade group and to the coupling between modes as the frequencies 
are close together. The mode shape of the fifth mode at zero revolutions per 
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second is similar to the mode shape of the sixth mode at 60 revolutions per 
second . 
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APPENDIX 


ANALYSIS 

r^f cccrr'iTT AL STlc::^!=SS ^r\Vtl A\n STATIC SOLUTION 
irPNT ICAl SU'^ST^^'lCTUCPS 

LISYINC pr PHAS*^ 1 

(INITIAI SUHStc^jCTuJ^E ANiLYSISI 


10 ^Icp^^n;),PH£$^ 

OIAG 2. a, 1^,14 

^PP risp 

TIMF 

S'^L i^.n 
Chkpnt v^s 
ALT pc AO, 40 

S'* A3 r,PI , KGGX/KGO/Vf N,LUSFT/VtNf NOCENL/V,N ,N0SlMOi 
ALTEO Ify 

LPLXf 
AlTfa 78 
LAH«^L LP'L X< 

ALTER «6 

F8S I oo,uon,pn/uccvi 

CHkPNT UOOVT 

OU^piJ'*'! K A!\,Pl *PAR,VECT! , P AR YFC T 3 , r» Ac^VF C T 3/ /C , N, - 1 /C * N » Ot 
OUT PUT 1 PAKVcCTA ♦f'ARVECTS , , •// C , N,0/C. N,OI 
ALTEa BY 145 
alter 151 15? 

ENDALTER 

CEND 

(CASE control DECKI 
BEGIN BULK 

-ll \rC-SSAPY STp TijPaL OATA PLUS THE SUHSTRUCTUR I NG 
'MATRIX PPEPA'^f^RS I 

ENPCATA 

END* % INCLUDE THIS CARO FOR 3^0/3T0 

rtVMENTS 

OISP APPROACH, 

ALL AriALYSIS SET DEGREES 0^ FREEOO^^ ShOULO BE INCLUDED ON ASET CAROS. 
PAPTITIOMING vectors which P-CVinf THE INFORMATION OF HOW THESE 
SUBSTR'JCTUPES ARE TIED TOGETHER* MUST PE INCLUDED IN BULK DATA CARDS. 
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F A^l/iLYSIS 

CIFF'^RENTI^I STIfft^SS MODAL AU-: STATIC SHLUTIOM 
!C*^ntical su^c'^cijrTuPC<; 

LISTI^G OP PHAS*" 1 ! 

(STATIC S'JBS^DiJCTUPP CniJPLPlG ANALYSISI 


DI cpYNM,PHACr; 

T I M r ^ -1 
APP OMAP 
OIAG 
p^r, 

PAP r.’A //r,ri,""^^/Vin,T'«uF = -l ^ 

INPIiTTl /KAA,PL, »./C*Nf-VCffM1 INPl 
PRE K AA=SAVr/PL = SAVF F 
LA‘Kl LCrf»911 

IKjPijTTi /P,, ,,/r,’^’,0/C*M*M 
M F P Gil f «f fKAAfFf/K GC T f 

AOO KGC.«KGGT/K^f 

FOUIV kT*kGG/TPmS1 

><FPGF* ♦PL. ♦ ♦ ♦E/PG’^/r .N,^U 

AOO PG.PGT/ptf 

eCHJlV PT,PG/TRIJC! 

pEPT LOnp9o*4,« 

PA«^N KGG.S^"CV,/KRFD.../C*Nt-lf 
PAP"^N PG, ,^Prv/PPEO. f./CtNtU 
SOLVE KCfO»DKeD/MLVT/C.N,lt 
maT»>PN ULVT. . . . //t 

vcoGEt UL VT, , , ♦ .SPr V/<ILVTT/C.N,1 $ 

MATPUN Ul yTT . . . .//T 

♦ WPITC MSfP'S *^AP- F^P PHAS=: S ^A^A pFCOVEOy. 

INPUTTI /,♦ . f /C fN.-3/C.NtU IM°1 

INPUTTI /♦,, ,/C,V,2/r,N,ii INPl 

OUTPUTS ♦ ♦♦ 

lA^EL LOCP?ft^ 

TNPdTTl /Q». ,,/F,Kt,0/CfNtll 
MATPPN Of ♦ •♦ //T 
PAPTM tjLVTTf,Q/,ULV,,/C.Nfli 
^^ATPPN UL V f f . . //^ 

f^llTrijTi C,iJLV..f//TlNPT 

PFPT LOOPP^.il 
OUTPljTl ♦ , ,♦ , y/G , ^'f-3! 

ENO! 

"=NJP 

(CASE CGNTurt n2r<» 

RFGIN PUl K 

MATRIX rr»rRATnps 

FNC rA"*A 

i INCLMOF THIS CAPO FDO IRV ?6J/^73 


OV*p APP'^OAGH# 

9FPGATING LOr’PS# 

ADOlTr*^NA( SINGL^ POINT CONS’^o A 1 NT5 AP H APPLIFP VIA MATF.lX PARTITION. 

partitionr'G vGfToF spcv mhst rf incluocO in bulk oata deck. 

The bulk data o^ck must include the omi capos fop the initiali- 

ZATIHN nc kGG ANO PG# 
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<;uBSTCl)rTiiRe ANALYSIS 

OlFPfRrjTlAl STIt^F^rsS M.^OAL -NO STATIC SPLUTICN 
ICfNrirM. SUHST?ijrT.)prs 
I ISTING pF PHAS‘D II I 

(STATIC DATA PPC'TVfPY AMO INITIAL PI^FtSrNTlAL STIFFNLSSI 


CIFFPYN'.PHAS'^B 
01 AG ?,e,I3,l4 
APp ni<;p 

SOL 13,0 

ChKPNT YfS 
TIME 7n 
ALTr^^ 3,“^ 

ALTf? 10,9^ 

iNPtjTTl /,, , ,/r ,M,-M 

nUTP'lTl, ,,, ,//C,':,-l/C ,N,1 1* 

PAD'M //r ,M,NnP/V,M,wACKr7t 
SAVE ^APKf 

PAPA'>i //C,N,NOP/V,'MH APF = 0 % 

JUMP L^npQCf 

LAPEL LOnp^^jf. 

PAFiM /AC,f , Ann/V,'.’,PLAOr /V ,N,HL A:if/C ,N ,1 i 
pCTPACM //C,N,1/C,N,GLAPF % 

INPUT^l /E,U1 V, , , /C 
AL^EP 103 

FILE <BLL=S AVf /MAA = «c AVF/PPL = SAVE f 
P AP AM //'“ ,N, SMB/V, N,MAO K/V, N,‘<A<1 K / : , N, U 
PPTPAOM //C f '‘1,0/r , V, vftt k $ 

CnNP 0 1FF3 ,‘.mpk! 

JUMP SKTOOFt 
LABEL D!Ff3f 

PAP AM //C ,N, ADC/V ,^',MAPK/V,N,MAF K/C,N, 100$ 

ALTf3 104 

SAVE rscrs^Tf 

105,105 
ALTFP 106,107 

SOU IV kOGG,KDNN/mpC^?/MGG,MNG/mdCF2$ 

alt-fp me, 10^ 

ALTFP 110, uo 

MCE? IISET,GM ,KOr,G,MGr,, ,/KCNN,MNN, ,1 

ALTER IIA,114 
ALTEP 116, U6 

SC El ilSFT,KDNN,MMN, ,/kOFP , kO^S , kOSS , MFF , , t 

ALTEP 117,117 


308 


n 







* 

It 







THKPNT KO^^S 1 
ALTCk wotiao 

AL^E=5 125 

»=0UIV t>L, PBl />^SCOSET/PS.P^<;/MSr OSE T/YS*Y0 S/OSCOS*^T/uGOV tUbOLV/ 

DSCOSET< 

ChkPNT P0L ,Pl3Sf YBStUBOOVi 

PAPAM //C,N,^PY/V,N,NDSK IP/C tNtO/'* 

OSMC? MPT,KAAfKrAA*KPStKOFS*K5SfKDS5tPL»PS,YSfU30V/KPLL,KBFS,KBSS, 
P^L .PpF.YbSf 1!9 00V/V»N fNOS*^ 1 P/V ,N ,OEPEATO/V* N, DSC JSETf 
<;avf N'Osk ip,p epeatp $ 

ChkPNT KbLL,KDf StKbSS tPPL tPeS*Yi S.UPOOV I 
lABEL SKfPOFl 

n’JTPijT] £t*,»//C*N,0/C,N,l$ 

PEPT LD0P99«A + 

nilTPuTi KbtL*^'AA,PPLf ,//C,N,0/C*N»11 

JUf^P FIN'ISF 

enoaltpp 

( 

f CH^CKFOINT nIrT^P^'APY ENT^P HF^P 

i 

C^ND 

(f ASP CP^.'TPnt rCCK ) 
per, I^J PMLK 
ENrOATA 

END* % INCLUDE THIS PDF I ?^0/'*TQ 




ApppQACH 01 SP. 

P P F T A P T . 

PEPFA^INO L^'OrS. 

T|4r n|cp»cprjTUL STIFFNESS *'AY CDNSI^tr^P AS r'»E\TlCAL 

FOP ALL ^UBSTMir Timf s 0^0710.-0 TMAT THF BIUND^PY »:PPPfTS A'i t NQT LA-^iiE. 
FOR saving CO»^PU^INr, TI'^F, THE CE'iTPk 0L AGE 0 I F r F r^T I AL STIPPNbSS 
P«ATCIX IS CHOSEN TO REPprsf^NT \LL. 

FOR OENESATP.n Th*: A-SET UI p :: w p NT I AL STIFFNP$S MATRIX, THE USER 
MAY riTHER rMOOSF USE SMPI Fip, SMP2, ThE lATEf- IS USCD IN 

THIS analysis, 

SOME nA*"A Sf'^F IN OHKPNT STa^P'^FNTS ARF SELECTIVELY DELETEDf ^UP 
NAS^PAN r:CE5 NOT ALLOW MJLTI-Rppl CHFCK-POINT TAPE, HPNCE CANNOT 
irrnvynOATP ALL ThF LAPDf siAf n^TA SETS# PPCGRAM INTL^RUPTKJN 
wnui 0 nCCilF IF yhF check-point TAPF had reached TO AN PNO. 
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« J'»S^PUC^i:=E ANALYSIS 

DTFPFPPNTIAl ^TtFr-NfSS MCC-L AN > ST^^IC SOLUTI'^fi 
I^ = N^1CAL ci|hSTPijrTtjPc$ 

LlSriMr, 0- PHt^r IV 

( DirpEPENTI AL ST!-^NESS STATIC ANO DYNA'MC ClUPLlNG ANALYSIS! 


I 

V 

? 

) 

i' 

) 



I 



10 OlP^DYNfOHASFA 
ni AG Z,^,13,1A,U 
App 01 sp 
T1M€ I )0 
SOL 13fO 

ALTtO 1 

PAPA^» //C»^»^OP/v,ri♦^=u^=-l % 

$ TMJ'^ USEO AS Di^A^STrP IN EC’JiV S"A"E^:NTS TO -OUtVALENCt '"ATA •’^LOCKS 

$ 

$ 0»^AP alter, sol 13tO PHASE IV, 

% RPPEA'^ING LOO'’. 

% C^'IT, SPC, ^‘PC A‘0 SUPC'ORT CARDS AC E PEP’^MTTEO wc;5P 

S 

4 USER MUS" USr SPOIfJT CARO TO JSE 0 ^ SPC AND Mor CA=DS 

4 USER !^’JS^ CREATE fJUU SQUA-E T AN: *^‘T ^^AToICcS «ITh ?mi cARCS IN 
4 '"UP bulk DA^A OECK 

4 

ALTPO 6*^1 
AL^ER A8,5o 

ALTEC 5<i.nO T SKIP static SOLUTION AND RORV'JLATION rc OIFFERCNTUL 

4 STI PCNESS NATC f X 

IVPUTTI 

INPUTTl /KC»c tVEL ,P^L**/C,N,*R/r ,N,0? SKIP « BLOCKS AND READ I^pT* 
pile CPP^SAVf /‘••BB = <AV?/PeL^$AVE 4 
INPUTTl /,,*»/C,N*-34 I\PT. 

JU'^P L^C^99 4 

LA«tL LOOP9^4 ANO ''l="l"I'Nt <=TIRPNESS A\0 *'‘A$S 'MATRICES 

IMPUTTl /Ef ♦ , ,/C ,M»04 c E LO P AP t j r HW IN G V EC T OR S FRQm I\PT. 

M^RGEf ••iKOfitff £f /KPAATli 

AOO KT *KntATl/KTT4 

EOUIV KTT,KT/TEM|Cf 

MERGE, , , ,MBe,E,/wiA'^l4 

APO mt,^AA"1/vTT4 

EOUIV MTT,MT/TFV^E4 

ME^GE, ,PBL, , , ,f /RBT/C ,N,4l4 

Apn PB,PPT/PT4 

POUtV PT,PR/TCuc|. 

0C9T LP0P99,«R4 TH^’Al 'MATRICES N Iw Fps'^cp k 

FOUIV ,kONN/*> FTF3 / HT,M^4s/M3r t> f 
CUNP LRL;0,^pCF? 1 
Mrci JSC^,PG/G*M 

,GM ,kT ,MT ♦ ,/<'^NN,*^VN, ,4 

alter 

sr=l •JSC’- ,F 0 MN,MNN, ,/< 0 f '<‘St‘^rP,,4 

alter 122,132 

S^p I IS^T tKPPP , , ,/G0,KOAA ,»NUL:. 1 .N ‘It ?,f»ULL ’ , , , , ,4 
AITPC 125 

PPMG2 KOAA/LLL,lJLL 4 
Chkpnt ULL 4 
EOUIV PR,PLM/NC$FT ♦ 
pmkpnt PL'^ f 

rONO PHA4 LI,^CSc^ 4 

SSG? USE1,GM,yS,‘‘TH^S,N0, ,p=^ / ,P"*3, p$ R, RIBI 
r.HK4NT Ppa,PSH,FL'’ 4 
LABEL PHA4LI 4 

SSni LLL ,Ul L •^DAA ,''LR,NULL: ,hULL ? ,N J L L 1 , >0 E*/ U Bl V , JB. V • R **UL V , ^ U' V/V,\, 

0^1T/V,Y, IP'^S.-'/r ,N,1 /V,N,EpSI 4 
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SAVP BPSl $ 

ruKT^.T URLVf t)0<^V ,?ULV , R'JQV ! 
ccf'P PH4AL?. es f 

W4TCC3 GPL ♦USPT♦'^TL, = ^^JLy//Cf^fL ! 

MfiTGpc r,9L , 1 »SET , ,p;<ijnv//C ,N, \ <■ 
l&pfL PHA4L2 1 

S?P1 USPT^P^i'Jbl'^ tU^OOV t/C*Nfl/C » 

fUPKLO « 

PURPE PBG^»f0^r,/TPJE ^ 

OAOAM //C »N* "lOP /V .M » = 0 ! 

iNPyTTl / * , , ,/C ,‘'.-3 t 

ouTPiiTi, , ,, ,//r •N*-i/c»p,,2 f pp: 

JUMP LOrP’^T f 
LAS'LL LT'^PO'^ t 

0\^tH //Cf N* ^US/V fN f S'H/C *\f 1 ! 

prtparm //:*\',d/c*m»S'J" f 
U:puTti /XV, , ,,/r ,^,0 1 

OiOTN UBT-Vf •XX/ , V, * /C f f I t 
matt CN UBV • f // ^ 

OUTPUT! XX,U6V,,,//C *‘1.0/C,N,? 1 lf^P2 
REPT L0QP97,4 ! 

OUTPUT!, ,,,,//r,^,-3/c,y,2 % P)P2 
AtT£W !‘»8,!45 
MA-^PCM PwtG, , , ,/ / $ 

OUTPUT ! , , ,, , //C • \,-l/C *N tl $INPl 

INPyTTl /,,, ,/<*,Nt-?i 

nuTPyTi L , , , , / /r. v,o/c ,^’ f I s 

PAC/v //r ,^,,NOP/V♦* .^LArp=0 ' 

Jlf'^P L00PP8 ♦ 

LABfL L'^OP98i PREPARE iNC PL-lc 3A®TITI^N=D = I C-pN VE C T S 0\ TiPE 

OA«A« //C ,* f 400/ V ,^I♦9U * ^E/ Vf ^, AT r/C f \ f 11 

POTPAP^ //'•,NtO/'**N,Bl.AOE 1 

INPUTTl /0,,,./C,M,0 f 

MATPPN 0. * f t //» 

PART^j pHlG♦fO/*PHl•,/^,^,l$ 

maTPPN pui , , , ,//l 

OUTPUT! PMI , , , , //C *NfO/C ,N , 1? 

REPT lO?P90f4 ! 

OUTPUT!,,,,, //r ,»i,-3/r ,N,i % 

ENDALtep 

C^NO 

(C A«^E C'^^TPOL “Frx \ 

B^GIN «ULx 

INCLUJ^ MATRIX f)0'“QATOPS 
INGL'ID? PPEUCO&TSuCTus ? CA-^A 
EN'^CAta 

FMO* $ I^CLU^E TUJ5 rAC^ 3'-J/3*’) 


-CC*^^E^iT5 

APOboich risp. 
pEPEATpiG L'^^P^. 

similar to tmas: II, BUT U? T ALT?c OIGIO PQ^^AT 13 lU PRrE^ TO 
EASE ThR |vPLr'•'‘^^ ATITN APPIXIVG SECPNO nPOE« 0-^IT ANC CP^ST®A!^4r5 
FOR '^|^,r^ *VALU^ TxTRifTi.vi, 

NO SECONP is APPlIEj in the PRcsb^t a\4LYSIS. 
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•*. S I ■’"• ■ .V *LV-! * 

1 rrrr ir ii <^-)* ^ r - i- = s 

LjQ-t.,o rr r--.c= V 

( A 5 =r V- - Y I 


9 

If 

•» 

t 

» 

f 

i 

> 





tr' =IVE! '.'•"('AL -.IccrLr - t . 

T I Mr 7 c 

APO !} I SP 
C I A C 2 • R • 1 3 , 1 ^ 

SOL !3,0 

20* l>5 

P&J3 AM / / : ♦ ‘i^VCP/V,:! J = = -l f 
I^PLT^l / , , , ,/r /C ,V*1 i 

jlivp I 

LA*’ 51 I A % 

ADO ii^ri r.A*PLU5iooo/T‘^Mp v 

EOUIV pjnic a/^?U5 5 

VATPPN I \niCA, , , ,// f 

OorrijT? 1»401 C A, , , , //^ /r *Y *n: - r^AC ! 

jmphtti /"pv, , , , / r ; /C •'•* 1 ! 

PIL*' ’J^nOVxS A Ve/Y^SrS *-V=/G^ = sA ;{ / M= ^ !7-/0’-'=S AV = /"‘hS5S S = SAV“ / 

<P^S-5i v«= I 

snc 1 USE Y , ,tj«v*'i‘onv ♦ YA$ , r,c •0'^ .r f** - p S*YH ss ♦ /uoov , * 

0^O/v*N ,N0E*^ I P/r *‘ *0$! I 
rM*(PM U?^»V*Qr<0 ? 

SDC2 CASeCC*CSYM,».rT,ru lL,r,PTY,i:pT, f^HPri, , 

^PGt’JnOv.^'ST *OU?OV! .0P^?l ,06Pbl*PU?Gvl/C*U*ubl % 

nfp nOPGl *nuPOVl*rcsp! .(TFrm, .//VtN.CAP'^NO s 

9fPT L^'^P9P*4 ♦ 

INPlfTTl /LAf^A*** •/C.^.-?f 
JUMP irOP^>QS 
tA^n LOnpeo^ 

A CO I\0ICA*PIU5 1000 /T= mo: f 

FOUIV Tf MPl * INC!CA/T?IJ5 » 

MATpPN INOKA,* • • // f 

0U^P'J"3 Tf.OTCA, ,*,//'• ,^',!/r, v,Ul= Vi-) ( 

tNPUTTl /OMtA,,*,/C*N*0i 
AlTt;* :C5 
pepY tr0P0P,4f 
cnOALTEP 

♦ C^•Er*<cnI^T '"ir'^l^SnCY =\TCS h?:e 

I 

f CASE C^NTsri, ) 

«^SGI\ P'JLK 
FNCTATA 

6Nf^^ S !NCLU^^ ^»-!S CA-D pv ifO/3'^0 


-•**Cr‘^vr'*TS- — 

40pp.-. ACM risp» 

QCOPA^IMr, LCCP', 

^M»S 0?».2'A*^PS A HijCr -c ST5FSS OUTPUT* A 0uTPuT3 

T- or ij^c- T- Micir ymp $’-“SS -Il^f niMCcwlSE IT WOUIO NOT 5 ?A$V 
TO pl*>y jt r* ^ '•AlCC^P pL"'^E 9 * 
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tHurn um^m ni l 




. ,J 5i i • t*' -t. 
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*■ ‘ ' 'iff ; 


nastran substructure analysis 

DYNAMIC SOLUTION NIThOUT 0 IFFE»PNTIAI STIFCNESS 
CMAP oROGBAM to CCMRInE TAPFS 

INPUT TAPF INPT CONTAINS STIFFNFSS, LOAD AND PAOTITink) MATRICES 
(PHASF I OUTPUT) 

INPUT fARE INPl CONTAINS MASS MATPIY 
IPHASE III OUTPUT) 

OUTPUT TAPE INP 2 



to TAPESfTv.rCNE ♦ 

TIME 2 
APP OMA» 

OIAG 2 , 8 . 13 , lA 
BEGIN f 

INPUTT I /KAA,PL,El,E2,E3/C,N,-3» 
INPUTTl /EA.E5, , , /C,N,0» 

INPUTTl /,,,,'C,N,-3/C,N,li 
INPUTTl /MAA,,,,/C,N,t'>/C.N, 11 
OUTPUTl, ,,, ,//C, N.-l/C ,N,2* 

OUTPUTl KAA,MAA,Ft,F2,E3//C,N,0/C,N,?t 
OUTPUTl FA,es,,,//C,N,0/C,N,2« 

OUTPUTl, , ,, ,//C, N.-3/C ,N,2I 
INPUTTl /,,,,/C,N.-3S 
END ♦ 

CENO 


ID BSVIRP,PHASE2 
time 95 
APP OISP 

OIAG 1,2,5,13,15,16 
SOL 3,0 
ALTER 1 
\ 

\ DMAP alter, sni 3,0 PHASE II. 

I REPEATING LOOP. 

» OMIT, SPC, MPC AO SUPPORT CAkDS ARE PEitMlTTEO HERE. 

\ 

PARAM //C,N,NOP/V,N,TRUE«-l» 

» TRUE USED AS PARAMETER IN EOUlV STATEMENTS TU EgUlWALENCE DATA RLJCNS 
PARAM //C,N,N0P/V,Y,ISTFSEE»-1 t 

% ISTFSFf crNTROLS WHETHER P1CT0«IAL lATktX PRINTER USED FOR STIFFNESS 

* TJ see USE ’>ARAM ISTFSEE I CARO IN BULK DATA DECK 

* MUST RF VARIABLE AS USEj IN CONO STATEMENT 
PARAM //C,N,N0P/V,Y,MASSSFE«-1 I 

» MASSSFF CONTROLS VHETHER PICTORIAL MATRIX PRINTER USED FOR MASS 
» TJ SEE USE PARAM MASSSfE 1 CARO IN BULK DATA DECK 

t MUST RE A VARIAJLE AS UScO IN CONO STATEMENT 

alter 6,41 

INPUTTl /KAAl ,MAA1, ,, /C ,N,>3/C,N,l t INPl — TnO TAPES 
I INPUTTl /Z,KAAl,MAAl,,/C,N,-3/C,N,l i INPl — SIK TAPtS 
FILE KAA1«SAVE/MAAI«SAVE» 

CONO SNSEEI.ISTFSEE % 

SEEMAT XAAl,,,,// t PRINTS LOCATION OF NCN-/LPO t ERMS OE KAAl MATRIX 
LAREL SMSFEl t 


i 


I 
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CONO f'^ISPClfMASSS^F i 

SEEMAT P»<INTS LJCATILN Cf- TFf^MS OF HAAl MATf^I^ 

LABEL ‘IMSEE; % 

% PARAM //CtN,NCP/V,N .IPT#0 S SIX TAPE'S 
LABEL LOno'^qs 
i BEGIN I )OP 99 

I PAR\^ //Cf N,A[)n/Vf%f IPT/V*Nf IPT/Cf\# I % SIX TA^FS 

% PRTPAR'- //C fN.O/C ,N • I PT S SIX TAPf;S 

$ INPllTTl / F # # . »/C f N#-3/V , I , IPT % SIX TAPES 

INPUTTI /f^ft t#/C,‘4,'J/CfNf I i INPI — Two TAPES 

HATPRN C* f f I //i 

MERGE , # • ,KAAl tE , /<S;Tl 

ADO KGS#KOGT/KTi 

I KT AND MT ARE Cf'NSI >h-irO \S SCRATCH r)ATA 3LJCKS AND MUST NfiT 
$ REFERENCED OUTSIDE JF LOOPS9 

FQUI V KT.KGG/TRUE % 

MFRGEf ,MAAltF */^w ,TS 
ADD MGG,Mj;T/MTf 
Fgulv MT , MGS/ TPijF V 
CONO SMSFF2, ISTFSEE » 

SEEMAT KGGft#*// t PRINTS LOCATI./, ‘IF NQN-ZEEO TERMS OF X GG MATRIX 

label SMSFE2 S 

CDND MMSE62.MASSSEE S 

SF6MAT MSGt*,,// I PRINTS LICATICN JF N)N-ZEhC TfPMS OF MGG MATRIX 
LABEL MMSEE2 % 

REPT L0JPO9,4S 

% THE 4 IN REPT LnDP99#*tl I'DICATFS TifAT LCJQP99 IS GONE THROUGH b TIMES 
$ Tn CHANGE NU^HFR .iF IJE\TIC*.L SJ3STHUCT JR ES ANALV/FD FROM St 

% CHANGF this NUMBER TO DNc LESS THAN THE NUMBER CF IDENTICAL 

% SUBSTRUCTURrS 

t END LSno 99 
ALTER 

alter 105»106 
OUTPUT I L AMA, ,, ,//r 

I PARAM //C#NfNrP/VfN,lPl# J I SIX TAPES 

INPUTTl /Ql,02tft/CfNt-3/C# 4,1 t I\Pl — Twj TAPES 

LABEL U.0P98I 

I PARAM //C.N,A^f)/V,N,lPl/V* M.IPl/C fNa S SIX TAPCS 
I PRTPAPM //CtN,J/r,N#l»l i SIX TAPFS 
% INPUTTl /O, ,,,/C, N.-3,V,N,tPl $ S?X TAPFS 
INPUTTl /O,,, ,/C,N, )/C,N, I i INPI -- TAPFS 

% Q CnRPFSPHNOS TO F IN L0Q®99 
PARTN PHIG,,Q/,PHI4,,/C,N,U 
OUTPUT I PHIZ,, . ,//: ,N,0/CtN,Oi 
REPT L30P9B,4t 

SO»2 CASECC,CSTv,>pt,DIT,p jpxIN,SIL • , , 3SPDT ,L AMA , GG ,PHIG , ,/, 
OQGl,nPM|G,,,/C,N,RtIGS 
OFP OPHIG,OOGl,,, ,//V,N,CAR )NOI 
alter 109.112 
alter UA,U5 
ENDALTFR 
CENO 


Calculation Times for Substructuring Analysis of 3 Lashed 
80-cm (31- in.) Steam Turbine Blades 





































TABLE 1 (Continued) 
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Figure 2 • Substructure Runs for Static or Dynamic 
(Natural Frequency) Analysis, with Differential Stiffness, 
of Identical Substructures. 







